Magnetic fields in circumstellar disks: The potential of Zeeman observations

Context. Recent high angular resolution polarimetric continuum observations of circumstellar disks provide new insights into their magnetic field. However, direct constraints are limited to the plane of sky component of the magnetic field. Observations of Zeeman split spectral lines are a potential approach to enhance these insights by providing complementary information. Aims. We investigate which constraints for magnetic fields in circumstellar disks can be obtained from Zeeman observations of the $113~\mathrm{GHz}$ CN lines. Furthermore, we analyze the requirements to perform these observations and their dependence on selected quantities. Methods. We simulate the Zeeman splitting with the radiative transfer (RT) code POLARIS (Reissl et al. 2016) extended by our Zeeman splitting RT extension ZRAD (Brauer et al. 2017), which is based on the line RT code Mol3D (Ober et al. 2015). Results. We find that Zeeman observations of the $113~\mathrm{GHz}$ CN lines provide significant insights into the magnetic field of circumstellar disks. However, with the capabilities of recent and upcoming instrument/observatories, even spatially unresolved observations would be challenging. Nevertheless, these observations are feasible for the most massive disks with a strong magnetic field and high abundance of CN/H. The most restrictive quantity is the magnetic field strength, which should be at least in the order of $\sim1~\mathrm{mG}$. In addition, the inclination of the disk should be around $60\deg$ to preserve the ability to derive the line-of-sight (LOS) magnetic field strength and to obtain a sufficiently high circularly polarized flux.Comment: 15 pages, 14 figure

Radiation Magnetohydrodynamics In Global Simulations Of Protoplanetary Disks

Our aim is to study the thermal and dynamical evolution of protoplanetary disks in global simulations, including the physics of radiation transfer and magneto-hydrodynamic (MHD) turbulence caused by the magneto-rotational instability. We develop a radiative transfer method based on the flux-limited diffusion approximation that includes frequency dependent irradiation by the central star. This hybrid scheme is implemented in the PLUTO code. The focus of our implementation is on the performance of the radiative transfer method. Using an optimized Jacobi preconditioned BiCGSTAB solver, the radiative module is three times faster than the MHD step for the disk setup we consider. We obtain weak scaling efficiencies of 70% up to 1024 cores. We present the first global 3D radiation MHD simulations of a stratified protoplanetary disk. The disk model parameters are chosen to approximate those of the system AS 209 in the star-forming region Ophiuchus. Starting the simulation from a disk in radiative and hydrostatic equilibrium, the magnetorotational instability quickly causes MHD turbulence and heating in the disk. For the disk parameters we use, turbulent dissipation heats the disk midplane and raises the temperature by about 15% compared to passive disk models. A roughly flat vertical temperature profile establishes in the disk optically thick region close to the midplane. We reproduce the vertical temperature profile with a viscous disk models for which the stress tensor vertical profile is flat in the bulk of the disk and vanishes in the disk corona. The present paper demonstrates for the first time that global radiation MHD simulations of turbulent protoplanetary disks are feasible with current computational facilities. This opens up the windows to a wide range of studies of the dynamics of protoplanetary disks inner parts, for which there are significant observational constraints.Comment: Accepted to A&

Radiation hydrodynamical models of the inner rim in protoplanetary disks

Many stars host planets orbiting within a few astronomical units (AU). The occurrence rate and distributions of masses and orbits vary greatly with the host stars mass. These close planets origins are a mystery that motivates investigating protoplanetary disks central regions. A key factor governing the conditions near the star is the silicate sublimation front, which largely determines where the starlight is absorbed, and which is often called the inner rim. We present the first radiation hydrodynamical modeling of the sublimation front in the disks around the young intermediate-mass stars called Herbig Ae stars. The models are axisymmetric, and include starlight heating, silicate grains sublimating and condensing to equilibrium at the local, time-dependent temperature and density, and accretion stresses parametrizing the results of MHD magneto-rotational turbulence models. The results compare well with radiation hydrostatic solutions, and prove to be dynamically stable. Passing the model disks into Monte Carlo radiative transfer calculations, we show that the models satisfy observational constraints on the inner rims location. A small optically-thin halo of hot dust naturally arises between the inner rim and the star. The inner rim has a substantial radial extent, corresponding to several disk scale heights. While the fronts overall position varies with the stellar luminosity, its radial extent depends on the mass accretion rate. A pressure maximum develops near the location of thermal ionization at temperatures about 1000 K. The pressure maximum is capable of halting solid pebbles radial drift and concentrating them in a zone where temperatures are sufficiently high for annealing to form crystalline silicates.Comment: accepted for Ap

Frontier Army Life Revealed by Charles King, 1844-1933

The purport of this study is to picture army life of the old West as revealed by the works of Charles King, author and military man of the late nineteenth and early twentieth centuries. I have emphasized the period of frontier policing by the army during the settling of the West, more commonly referred to as the Indian Wars of the 1870\u27s.https://scholars.fhsu.edu/fort_hays_studies_series/1023/thumbnail.jp

Gaps, Rings, and Non-Axisymmetric Structures in Protoplanetary Disks - From Simulations to ALMA Observations

Recent observations by the Atacama Large Millimeter/submillimeter Array (ALMA) of disks around young stars revealed distinct asymmetries in the dust continuum emission. In this work we want to study axisymmetric and non-axisymmetric structures, evocated by the magneto-rotational instability in the outer regions of protoplanetary disks. We combine the results of state-of-the-art numerical simulations with post-processing radiative transfer (RT) to generate synthetic maps and predictions for ALMA. We performed non-ideal global 3D MHD stratified simulations of the dead-zone outer edge using the FARGO MHD code PLUTO. The stellar and disk parameters are taken from a parameterized disk model applied for fitting high-angular resolution multi-wavelength observations of circumstellar disks. The 2D temperature and density profiles are calculated consistently from a given surface density profile and Monte-Carlo radiative transfer. The 2D Ohmic resistivity profile is calculated using a dust chemistry model. The magnetic field is a vertical net flux field. The resulting dust reemission provides the basis for the simulation of observations with ALMA. The fiducial model develops a large gap followed by a jump in surface density located at the dead-zone outer edge. The jump in density and pressure is strong enough to stop the radial drift of particles. In addition, we observe the generation of vortices by the Rossby wave instability (RWI) at the jumps location close to 60 AU. The vortices are steadily generated and destroyed at a cycle of 40 local orbits. The RT results and simulated ALMA observations predict the feasibility to observe such large scale structures appearing in magnetized disks without having a planet.Comment: Language update, added comments, added citations, in press. (A&A

3D MHD Simulations of Planet Migration in Turbulent Stratified Disks

We performed 3D MHD simulations of planet migration in stratified disks using the Godunov code PLUTO, where the disk is turbulent due to the magnetorotational instability. We study the migration for planets with different planet-star mass ratios $q=M_{p}/M_{s}$. In agreement with previous studies, for the low-mass planet cases ($q=5\times10^{-6}$ and $10^{-5}$), migration is dominated by random fluctuations in the torque. For a Jupiter-mass planet $(q=M_{p}/M_{s}=10^{-3}$ for $M_{s}=1M_{\odot})$, we find a reduction of the magnetic stress inside the orbit of the planet and around the gap region. After an initial stage where the torque on the planet is positive, it reverses and we recover migration rates similar to those found in disks where the turbulent viscosity is modelled by an $\alpha$ viscosity. For the intermediate-mass planets ($q=5\times10^{-5}, 10^{-4}$ and $2\times10^{-4}$) we find a new and so far unexpected behavior. In some cases they experience sustained and systematic outwards migration for the entire duration of the simulation. For this case, the horseshoe region is resolved and torques coming from the corotation region can remain unsaturated due to the stresses in the disk. These stresses are generated directly by the magnetic field. The magnitude of the horseshoe drag can overcome the negative Lindblad contribution when the local surface density profile is flat or increasing outwards, which we see in certain locations in our simulations due to the presence of a zonal flow. The intermediate-mass planet is migrating radially outwards in locations where there is a positive gradient of a pressure bump (zonal flow).Comment: Accepted for publication in Ap

Intrinsic polarisation of elongated porous dust grains

ALMA observations revealed recently polarised radiation of several protoplanetary disks in the (sub-)millimetre wavelength range. Besides self-scattering of large particles, thermal emission by elongated grains is a potential source for the detected polarisation signal. We calculate the wavelength dependent absorption and intrinsic polarisation of spheroidally shaped, micrometre and sub-millimetre sized dust grains using the discrete dipole approximation. In particular, we analyse the impact of dust grain porosity which appears to be present in disks when small grains coagulate to form larger aggregates. For the first time our results show that (a) the intrinsic polarisation decreases for increasing grain porosity and (b) the polarisation orientation flips by 90 degree for certain ratios of wavelength to grain size. We present a new method to constrain grain porosity and the grain size in protoplanetary disks using multi-wavelength polarisation observations in the far-infrared to millimetre wavelengths. Finally, we find that moderate grain porosities ($\mathcal{P}\lesssim0.7$) potentially explain the observed polarisation fraction in the system HD 142527 while highly porous grains ($\mathcal{P}>0.7$) fail unless the grain's axis ratio is extraordinarily large.Comment: 10 pages, 10 figure